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Abstract 

Solutions to flat space Friedmann-Robertson- Walker cosmologies in Brans- 
Dicke theory with a cosmological constant are investigated. The matter is 
modelled as a 7-law perfect fluid. The field equations are reduced from fourth 
order to second order through a change of variables, and the resulting two- 
dimensional system is analyzed using dynamical system theory. When the 
Brans-Dicke coupling constant is positive (u > 0), all initially expanding 
models approach exponential expansion at late times, regardless of the type 
of matter present. If u < 0, then a wide variety of qualitatively distinct models 
are present, including nonsingular "bounce" universes, "vacillating" universes 
and, in the special case of uj = — 1, models which approach stable Minkowski 
spacetime with an exponentially increasing scalar field at late times. Since 
power-law solutions do not exist, none of the models appear to offer any 
advantage over the standard deSitter solution of general relativity in achieving 
a graceful exit from inflation. 
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I. INTRODUCTION 



A recent renewal of interest in Brans-Dicke (BD) theory Pj] can be traced to the discovery 
by La and Steinhardt that the use of BD theory in place of general relativity can ameliorate 
the exit problem of inflationary cosmology @ . This is possible because the interaction of the 
BD scalar field with the metric slows the expansion from exponential to power-law. Although 
the original "extended inflation" scenario appears to have been ruled out by observational 
constraints [BJ, models which survive include those based on more general scalar-tensor 
theories, such as "hyperextended inflation" [|J, and hybrid models which include both a first- 
order phase transition and a period of slow- roll The renewed interest in scalar-tensor 
gravitation has led to several recent investigations into the generation of exact solutions for 
cosmology in such theories as well as to some qualitative studies of the models which 
result [|7|-|12|. It has also been pointed out recently that an inflationary era may result directly 
from the dynamics of the scalar field, without any potential or cosmological constant being 
necessary [|I3| . 

In this paper we are concerned with the behavior of homogeneous and isotropic cosmo- 
logical models in Brans-Dicke gravity, with the addition of a positive cosmological constant. 
This differs from standard extended inflationary scenarios in that the vacuum energy is de- 
coupled from the scalar field. The goal is simply to analyze the cosmological models which 
such a theory gives rise to, with an emphasis on the question of whether a viable inflationary 
model might exist. Previously, similar analyses have been performed for the case A ^ with 
no matter present ||10|| , and for the case A = with additional matter present |IT| , |12|| . The 



treatment will closely parallel that given in Ref. |12| , and the interested reader is referred 
to that paper for more detail. 

In Sec. II, it is shown that the field equations for this theory can be reduced to a two- 
dimensional dynamical system in the case of flat space. In Sec. Ill, the equilibrium points 
are found and the corresponding solutions are discussed. In Sec. IV we summarize the 
results. 



II. THE FIELD EQUATIONS 

In this section the field equations are reduced to a planar dynamical system through 
a change of variables. For notation and conventions, the reader is referred to the parallel 



treatment given in [ 12| . Generalizing the action for Brans-Dicke theory to include a nonzero 



cosmological constant, it may be written as 

Sbd = J d A x^ (-<f>[R - 2A] + + 167r£ m ^ . (2.1) 

Taking A and u constant, and varying this action with respect to the metric and the scalar 
field, one finds that the nontrivial components of the field equations in a homogeneous and 
isotropic (FRW) spacetime are 

/ • v 2 

(2.2) 
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where a(£) is the cosmic scale factor. Assuming a perfect fluid form for the stress-energy 
tensor, i.e., T^ u = diag(p, p,p,p), the usual conservation equation is also satisfied (the 
zeroth component of T^ h '. u = 0): 

p = -3^ + p). (2.4) 



Assuming only that p > and <p > 0, inspection of Eq. ( |2.2| ) reveals that we must have 
uj > —3/2 in order to satisfy the field equations for all values of k and A. Furthermore, we 
see from the form of the action in Eq. ( |2.1| ) that the integrity of the theory is lost if u — 0. 
We therefore take uj > —3/2 and uj ^ in what follows. 

Now we take k — 0, and transform the fourth-order system specified by Eqs. ( |2.2| - [2.4| ) 
into a pair of coupled second-order equations in which, however, only first derivatives of the 
new variables appear. It will be sufficient for our purposes to consider only the models with 
k = 0, as any candidate for a viable inflationary model must at the very least solve the 
flatness problem. First, define the new variables 

HM). <«> 

S = A^, (2.6) 

where dots represent derivatives with respect to time, and A = (2u + 3)/12. Next, 
parametrize the equation of state by writing p = (7 — l)p, where, for example, 

0, false- vacuum energy; 

■ 1, pressureless dust; , , 

4/3, radiation; 

2, "stiff" matter. 



Finally, take k = and rewrite Eqs. ( |2.2| - |2.4| ) in terms of these new variables. Then straight- 
forward differentiation and resubstitution lead to the equivalent field equations 

t = 2 (1 - 3 7 /4) - ^(1 - 3 7 /2) - 30S - ^(1 - 3 T /2), (2.8) 

Eqs. (|2.8|) and (|2.9|) constitute a planar dynamical system in the variables and S, and are 
the desired results of this section. 



III. THE EQUILIBRIUM POINTS 

The equilibrium points of the dynamical system are obtained by setting and £ equal 
to zero in Eqs. (|2.8|) and (|2.9|), and then solving the resulting equations for and E. As 
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each of these equations is a second order polynomial, we expect in general four equilibrium 
points. When = S = 0, this implies that a/a = — H/2A and 0/0 = T,/A are both 
constants, so that 



a(t) 
0(t) 



a exp 



^oexp 



0- 

Et 
~A 



2A 



(3.1) 
(3.2) 



Therefore a fixed point in the O-S system represents deSitter spacetime, with the addition 
of an exponentially varying scalar field. In the special case = S/2A, the solution is not 
deSitter but rather Minkowski spacetime with a scalar field. Note also that only equilibrium 
points at finite values have been considered. The global picture, including the behavior 
of and S at infinity, may be obtained by various compactification methods [13 j. Such 
an analysis is not necessary for our purposes, however, as we are primarily interested in 
whether viable inflationary models exist in this theory. In particular, one sees by inspection 
of Eqs. ( |2.8| ) and Q2.9|) that the origin of the 0-S plane can never be an equilibrium point 
when A ^ 0. This immediately rules out the possiblity of a stable power-law solution, since 
such a solution would appear in that plane as line of constant slope, and would asymptotically 
approach equilibrium at the origin. 

Although solution curves span the entire 0-S plane, the requirement p > 0, where p is 
the energy density of the perfect fluid matter, eliminates some regions on physical grounds. 
It follows from Eq. ( |2.2| ) with k = that (assuming > 0) any point (0o, So) satisfying 
©o 2 > Sq/A + A/3 will lie in p > 0, whereas points satisfying O 2 < SqA4 + A/3 lie in 
p < and thus do not represent physical solutions. Now let us proceed with the analysis. 
Although we restrict ourselves to the consideration of models with A > 0, it is clear that 
the techniques can easily be extended to models with a negative cosmological constant. 

One pair of equilibrium points are always present regardless of the value of 7; they satisfy 
the field equations for p = and are thus vacuum solutions. These points are 



(©o, s ; 



1.2 



± 



A(2cj + 3) 
2(3lu + 4) 



.1/2 



:m/6) 



(3.3) 



and they represent the solutions 

a(t) = aoexp |±(cj + 1) 

0(i) = 0o exp < ± 



2A 



(2w + 3)(3a; + 4) 



1/2 



2A 



(2cu + 3)(3a; + 4) 



1/2 



(3.4) 
(3.5) 



These solutions have previously been noted in the literature |7|JT0[1, and are attractors for 
most, but not all, of the initially expanding models in this theory (cf. discussion below). 
Note that if uj = —1, then there are solutions where the geometry is Minkowski and the 
scalar field either grows or shrinks exponentially. Under the field redefinition — > e~*, these 
can be identified as the static "linear dilaton" solutions of string cosmology [ITS . 



The location and stability of the remaining equilibrium points depends upon the value 
of 7, as well as upon uj and A. Hence it is convenient to classify the models, and to discuss 
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the overall character and stability of the solutions, according to the equation of state of the 
perfect fluid. The additional equilibrium points follows: 




(3.6) 
(3.7) 

(3.8) 
(3.9) 



The details of the stability analysis are given in Appendix A. Table 1 of that appendix 
lists the eigenvalues of the equilibrium points, and Table 2 explicity states the existence and 
stability of each equilibrium point as a function of uj. Note that in some cases, the solutions 
represented by a given equilibrium point require negative energy density, and are therefore 
physically uninteresting. In such cases we have simply written "p < 0" in the tables. 

The overall character of the solutions may be further examined by numerically integrat- 
ing the solutions 0(i) and S(t) for a variety of initial conditions with each qualitatively 
distinct set of parameters (A, a;, 7). Figures 1-4 show the results of this procedure, where 
we have selected A = 3 arbitrarily. The shaded regions in each case require p < 0, and so 
are disallowed physically. Curves to the right of the line da/dt = represent expanding 
universes, and those to its left represent contracting universes. Note that if — 1 < uj < 0, as 
for example in Fig. lb, then there exist nonsingular "bounce" models which pass smoothly 
from contraction to expansion, and "vacillating" models which pass from expansion to con- 
traction to reexpansion, or the time reversal of this behavior. There also exist extreme cases 
of "vacillation" ; some models can vacillate several times before settling down to exponential 
contraction at late times, as shown in Fig. 4d. 



IV. SUMMARY 

We have shown that in Brans-Dicke theory with a positive cosmological constant, a wide 
range of flat-space models exist, including some with no analogues in general relativity. In 
general, these models are parametrized by the initial conditions of the scalar field, the value 
of the BD coupling constant, an initial expansion rate, an equation of state for the matter 
and the value of the cosmological constant. The first two parameters are, of course, absent 
in general relativity. Most, but not all, of the expanding models asymptotically approach 
vacuum deSitter spacetime at late times. Power-law expansion is not possible when A is 
nonzero. 

If uj > —1, there exist two finite- valued equilibrium points of physical interest, repre- 
senting the vacuum solutions ( |3.4| ) and (|3~5|). One of these is stable and corresponds to 
expanding deSitter spacetime, and all initially expanding models approach this solution 
asymptotically. As uj — > 00, the solution becomes identical to the deSitter universe where 
a(t) ~ e At//3 and (p(t) = constant, in accordance with the well-established correspondence 
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between GR and BD theory in this limit. When lo > 0, all contracting models contract to 
a singularity; however, if — 1 < uj < 0, then nonsingular "bounce" models are also possible, 
and these may or may not approach deSitter spacetime at late times, depending upon the 



initial conditions of the model. This behavior has been noted by other authors [[U|. Also 



there are "vacillating" models, which expand from a big bang, slow down, and recontract 
before continuing their expansion and approaching deSitter spacetime. The time-reversal of 
this behavior also exists. 

If u — — 1, there exist "static-exponential" solutions, where the geometry is Minkowski 
spacetime while the scalar field changes exponentially with time. The solution with increas- 
ing scalar is found to be stable, and all models which expand from a big bang approach it 
at late times, regardless of the type of matter present. The stability of these models may 
be explained by the fact that in the Newtonian limit, ~ G^ 1 [fL7[]. Hence the exponential 
increase in the scalar field corresponds to an exponential weakening of the gravitational 
interactions, ensuring that the universe does not recollapse regardless of its matter content. 

If —3/2 < to < — 1, then the behavior of the models depends upon the type of matter 
present, and we can distinguish two classes of behavior, (i) In the cases of false-vacuum 
energy and pressureless dust, only the equilibrium points representing the vacuum solutions 
( ^.4| ) and ( |3.5| ) are present in the regime p > 0, and the contracting solution is stable. 
All models which start from a big bang eventually contract to a singularity; this collapse 
will become asymptotically exponential if —4/3 < u < — 1, or superexponential if —3/2 < 
uj < —4/3. There also exist models which start with a finite rate of expansion and expand 
perpetually, (ii) In the cases of radiation and "stiff" matter, the dynamical system undergoes 
a qualitative change at a particular value of uj. This critical value is uj c = —5/4 for radiation, 
and uj c = —7/6 for "stiff" matter. If uj c < uj < —1, then only the equilibrium points 
representing ( |3T4] ) and ( |3.5| ) are present with p > 0, and the contracting solution is stable. 
If —4/3 < uj < uj c , then there are four equilibrium points in the regime p > 0, representing 
both vacuum and non-vacuum deSitter solutions; however, only the contracting non- vacuum 
solution is stable. In these solutions, the exponential growth of the scalar field exactly 
balances that of the energy density, so that the ratio p/<f) is constant. Thus the ordinary 
matter acts exactly like a cosmological constant, since it is this ratio which appears as the 
matter source term in the field equation (|2.2|) . All models will collapse to a singularity; 
the collapse will be asymptotically exponential for models starting from a big bang. If 
—3/2 < uj < —4/3, the unstable vacuum equilibrium points no longer exist; otherwise the 
behavior of the models remains unchanged. 

Although somewhat exotic, the models discussed here do not seem to have any hope of 
solving the graceful exit problem of inflationary cosmology. In models of extended inflation, 
the mediation of the scalar field slows the expansion from exponential to power-law, so 
that the Hubble parameter decreases with time and true-vacuum bubble nucleation may 
complete the inflation-ending phase transition. Here, however, the cosmological constant 
induces deSitter spacetime, and all of the problems of Guth's "old" inflation recur. In 
addition, one is faced with the question of the origin of the cosmological constant in these 
models. Although there is no a priori reason to exclude such a term from the field equations, 
the usual explanation of a field with a nonzero potential as the source of the vacuum energy 
is not available in this case, since such a potential term would be coupled to the scalar field. 
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APPENDIX A: 

Defining £W = — @ and £ < - 2 - ) = S — X , the dynamical system specified by (|2.8| 



and ( |2.9| ) may be written in the form 



^ = J? + ••• , (Al) 



dt 

where the Jacobian is 

_ / -3 7 6 + S /2A -6(7/2 - l)E /A + e /A 

Ue (l - 37/4) - 3S -(1 - 3 7 /2)E /A - 30 o 



(A2) 



In cases where the eigenvalues of the Jacobian all have nonvanishing real part, the fixed 
point is called hyperbolic and we can determine its stability from the signs of those real 
parts: if the real part of each of the eigenvalues is negative at a given equilibrium point, 
the solution is stable at that point; if the real part of each eigenvalue is positive, or if the 
real part of one eigenvalue is positive and that of the other is negative, then the solution is 
unstable at that point. Finally, if the real part of any of the eigenvalues is zero at a point, 
then the point is called nonhyperbolic and its stability in the neighborhood of that point 
cannot be determined by this method [14[ . 



Table 1 shows the eigenvalues of the Jacobian for each equilibrium point, and Table 2 
explicitly states the existence and stability of the equilibrium points as a function of uj. 
In cases where the solution represented by the equilibrium point requires negative energy 
density, we have simply written "p < 0". The points are labelled in accordance with the 
conventions in the text, i.e., for each value of 7, points 1 and 2 represent vacuum solutions, 
and points 3 and 4 represent non-vacuum solutions. 
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TABLES 



TABLE I. Eigenvalues of the Jacobian Matrix for BD Cosmology with A > 



Equilibrium Point Ai A2 



(7 


= 0)! 


-[2A/(2uj + 3)(3kj + 4)] 1 /2 


-[2A(3u; + 4)/(2lj + 3)] 1 / 2 


(7 


= 0) 2 


+ [2A(3u; + 4)/(2^ + 3)] 1 /2 


+ [2A/(2uj + 3)(3uj + 4)] 1 /2 


(7 


= 0)3 


p < 


p < 


(7 


= 0)4 


p < 


p < 


(7 


= l)i 


-[2A(3^ + 4)/(2cj + 3)] 1 /2 


-[2A(3^ + 4)/(2cj + 3)] 1 /2 


(7 


= 1)2 


+ [2A(3cu + 4)/(2l^ + 3)] 1 /2 


+ [2A(3c^ + 4)/(2^ + 3)] 1 /2 


(7 


= 1)3 


p < 


p < 


(7 


= 1)4 


p < 


p < 


(7 = 


= 4/3)! 


-[2A(3^ + 4)/(2cj + 3)] 1 / 2 


-(4u + 5)[2A/(2cj + 3)(3w + 4)] 1 / 2 


(7 = 


= 4/3)2 


+ (4lo + 5)[2A/(2w + 3)(3w + 4)] x /2 


+ [2A(3c^ + 4)/(2^ + 3)] 1 /2 


(7 = 


= 4/3)3 


C(-l + V64w + 81) a 


C(-l - V64u; + 81) a 


(7 = 


= 4/3)4 


C(+l + V64w + 81) a 


C(+l - V64u; + 81) a 


(7 


= 2)i 


-[2A(3u; + 4)/(2i i j + 3)] 1 / 2 


-(6w + 7)[2A/(2cj + 3)(3w + 4)] 1 / 2 
+ [2A(3u; + 4)/(2 W + 3)] 1 /2 


(7 


= 2)2 


+(6u + 7)[2A/(2w + 3)(3w + 4)] x /2 


(7 


= 2)3 


C(- V3(2w + 3) + y/102u + 121) a 


C(-v / 3(2o; + 3) - >/102cd + 121) a 


(7 


= 2)4 


C(+V3(2u; + 3) + y/W2u> + 121) a 


C(+v/3(2cj + 3) - V102W + 121) a 


a C= [A/2(2w + 3)] 1 /2 



TABLE II. Existence and Stability of the Equilibrium Points 



Equilibrium Point 


-3/2 < lo < -4/3 


-4/3 < u < uj c a 


UJ > LO c 


(7 = 0, l)i 


nonexistent 


N/A 


stable 


(7 = 0, 1)2 


nonexistent 


N/A 


unstable 


(7 = 0, 1)3 


p < 


N/A 


p < 


(7 = 0, 1)4 


p < 


N/A 


p < 


(7 = 4/3, 2)x 


nonexistent 


unstable 


stable 


(7 = 4/3, 2) 2 


nonexistent 


unstable 


unstable 


(7 = 4/3, 2)3 


stable 


stable 


p < 


(7 = 4/3, 2)4 


unstable 


unstable 


p < 


a u; c (7 = 0, 1) = -4/3 


; c c ( 7 = 4/3) = -5/4 ; w c ( 7 = 2) 


= -7/6 
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FIGURES 



FIG. 1. The evolution of solutions in the 0-S plane with 7 = (false- vacuum energy), where 
A = 3 has been selected for convenience, (a) 10 = 9/2, representative of uj > 0. Solutions to 
the right of the line da/dt = represent expanding universes; those to the left are contracting 
universes. The shaded region requires negative energy density and so is disallowed physically, (b) 
uj = —1/2, representative of — 1 < uj < 0. The line da/dt = has moved into the regime of positive 
energy density, so that some models pass smoothly from contraction to expansion, a general feature 
of models with uj < 0. (c) uj = — 1. Note the attractive nature of the "static-exponential" solution, 
(d) uj = —1.4, representative of —3/2 < uj < —4/3. All models collapse to a singularity regardless 
of initial conditions. 

FIG. 2. Some models with 7 = 1 (pressureless dust) and A = 3. (a) uj = 9/2, representative of 
uj > 0. (b) uj = —1/2, representative of — 1 < uj < 0. (c) uj = —1. (d) uj = —1.4, representative of 
-3/2 < u < -4/3. 

FIG. 3. Some models with 7 = 4/3 (radiation) and A = 3. (a) uj = —1/2, representative of 
-1 < uj < 0. (b) uj = -I. (c) uj = -1.2, representative of -5/4 < uj < -1. (d) u = -1.4, 
representative of —3/2 < uj < —5/4. 

FIG. 4. Some models with 7 = 2 ("stiff matter") and A = 3. (a) uj = —1/2, representative 
of — 1 < uj < 0. (b) uj = —1. (c) uj = —1.2, representative of —7/6 < uj < —1. The loca- 
tions of the equilibrium points are marked with dots for clarity, (d) uj = —1.4, representative of 
— 3/2 < uj < —7/6. Only one solution is shown to simplify the plot. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9409002v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9409002v2 
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